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Carbon chemistry in Galactic Bulge Planetary Nebulae 



o 

(N 

X) 
D 

CnI 
(N 

<: 
O 

^• 

Oh! 

O 



oo 
O 



L. Guzman-Ramirez/* A. A. Zijlstra/ R. Ni Chuimm/ 
K. Gesicki,2 E. Lagadec,^ T. J. Millar^ and Paul M. Woods ^ 

^Jodrell Bank Centre for Astrophysics, School of Physics & Astronomy, The University of Manchester, Manchester, M13 9PL, UK 

^ Centrum Astronomii UMK, ul.Gagarina 11, 87-100 Torun, Poland 

^European Southern Observatory, Karl-Schwarzschild-Str 2, 85748 Garching, Germany 

"^Astrophysics Research Centre, School of Mathematics and Physics, Queen 's University Belfast, Belfast BT7 INN, UK 



Accepted 2010 November . Received 2010 November 04 



ABSTRACT 

Galactic Bulge Planetary Nebulae show evidence of mixed chemistry with emission from both 
silicate dust and PAHs. This mixed chemistry is unlikely to be related to carbon dredge up, 
as third dredge-up is not expected to occur in the low mass Bulge stars. We show that the 
phenomenon is widespread, and is seen in 30 nebulae out of 40 of our sample, selected on 
the basis of their infrared flux. HST images and UVES spectra show that the mixed chemistry 
is not related to the presence of emission-line stars, as it is in the Galactic disk population. 
We also rule out interaction with the ISM as origin of the PAHs. Instead, a strong correlation 
is found with morphology, and the presence of a dense torus. A chemical model is presented 
which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase 
chemical reactions. The model predicts two layers, one at Ay ~ 1 .5 where small hydrocarbons 
form from reactions with C^, and one at Ay ~ 4, where larger chains (and by implication, 
PAHs) form from reactions with neutral, atomic carbon. These reactions take place in a mini- 
PDR. We conclude that the mixed chemistry phenomenon occurring in the Galactic Bulge 
Planetary Nebulae is best explained through hydrocarbon chemistry in an UV-irradiated, dense 
torus. 
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1 INTRODUCTION 

Planetary Nebulae (PNe) are remnants of the extreme mass loss ex- 
perienced by Asymptotic Giant Branch (AGB) stars, a phase that 
occurs during the late stages of evolution of low and intermediate- 
mass stars, between 1 and 8 Mg. The intense mass loss (from 10"^ 
to 10"'* Mo yr"') leads to the formation of a circumstellar envelope 
made of gas and dust (Kwok 2000). The dust strongly emits in 
the infrared. After the termination of the AGB, the star rapidly in- 
creases in temperature and photodissociates and photoionizes the 
expanding ejecta. 

During the AGB phase, a star may evolve from being oxygen- 
rich to being carbon-rich, a change reflected in the chemical com- 
position of the stellar wind. The change occurs when carbon pro- 
duced by He-burning is brought to the surface by processes in the 
stellar interior thereby increasing the C/O ratio until it exceeds 
unity and a carbon star is formed. This is the so-called third dredge- 
up, w hich occurs at the end of flash burning in t he He-shell I Herwis 
l2005l) . The process depends on stellar mass: IVassiliadis & Wood 
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( Il993h showed the third dredge up occurs if the core mass of the 
star M„ > 1.5 Mq. 

This predicts a clear distinction, with some (lower-mass) stars 
showing oxygen-rich and some (higher-mass) stars carbon-rich 
ejecta. In the molecular ejecta, the CO molecule locks away the 
less abundant element, leaving the remaining free O or C to drive 
the chemistry and dust formation. Oxygen-rich shells are charac- 
terized by silicate dust, while carbon-rich shells show Polycyclic 
Aromatic Hydrocarbon (PAH) emission bands and carbonaceous 
dust. This dichotomy is indeed by and large observed. 

The first evidence for mixed-chemistry in post-AGB objects 
was found when one PN (IRAS07027-7934) with strong PAH emis- 
sion bands at 6.2, 7.7, 8. 6, and 11.3 /^m, wa s found to also show 
a 1.6GHz OH maser line i jZijlstra et al.|l99l |). Observations made 
with ISO uncovered several further cases, where PAH emission in 
PNe occurred together with emission bands of silicates at 23.5, 
27.5, 33.8, and 10 fim, usually found in 0-rich shells (Waters et 
al. 1998a,b; Cohen et al., 1999, 2002). 

The mixed- chemistry phenome non has also been seen in the 
Red Rectangle jWaters et al]|l998l) and in PNe with late Wolf- 
Rayet type central stars jWaters et al.lll998l) , where it is attributed 
to the presence of old oxygen-rich material in a circumstellar disk, 
with the PAHs forming in a more recent carbon-rich outflow. Weak 
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PAH emission has however been seen in som e clearly oxygen-ric h 
objects, such as NGC 6302 and Roberts 22 jZiilstra et atlbOOlh . 
PAH bands have been detecte d in S-type AGB st a rs, vyith C/0 ratios 
close to, but below, unity ( Smolders et al ] l2O10l) .l Cohen & Barlowl 
( l2005h have shown that some PAH emission can occur at C/O ratios 
as lo w as 0.6. 

IPerea-Calderon et al.l ( l2009l) showed that the mixed-chemistry 
phenomenon is widespread among the PNe in the Galactic Bulge 
(GB). Their Spitzer observations show that the simultaneous pres- 
ence of oxygen and carbon-rich dust features is common, and is 
not restricted to objects with late/cool [WC] type stars. The tradi- 
tional explanation relating the mixed chemistry to a recent evolu- 
tion towards carbon-star is unlikely for the Bulge objects, as these 
old, low-mass stars are not expected to show third dredge-up, and 
therefore should not show enhanced C/O ratios. The few AGB 
carbon stars in the Bu lge do not originate from third dredge- up 
jAzzopardi et al.|[l988[) . Different explanations are needed. 

In this paper, we investigate two alternative scenarios for the 
presence of C-rich molecules: PAHs swept-up from the interstellar 
medium (ISM), and in-situ PAH formation via chemical pathways 
in oxygen-rich environments. This is performed through the anal- 
ysis of 40 PNe towards the Galactic Bulge for which Spitzer spec- 
troscopy is available. HST images are available for 22 of these. 
In Section 2 we present the observations. In Section 3 we present a 
correlation between the silicate emission and the presence of a torus 
in the PNe from HST images. In Section 4 we present a chemical 
model that forms long C-chains in an O-rich environment and in 
Section 5 we discuss these results and present our main conclu- 
sions. 



2 OBSERVATIONS 

iPerea-Calderon et al. I tog*) uncovered the mixed chemistry in 2 1 
Galactic Bulge PNe, out of a total sample of 26, using Spitzer Space 
Telescope spectra. We retrieved and re-reduced these 21 Spitzer 
spectr a. We also included the PN H2-20 from iGutenkunst et al.l 
( l2008h . and a further 1 8 PNe in the GB for which we found spec- 
troscopic data in the Spitzer archive, yielding a total sample of 40. 

All the spectra wer e taken with the Infrared Spectrograph 
(IRS) faouck et alj|2004) . The spectra cover a range between 5.2- 
37.2 by using the Short-Low (SL: 5.2-14.5 yum; 64<R<128), 
Short-High (SH: 9.9-19.6 yum; R~600) and Long-High (LH: 18.7- 
37.2 yum; R~600) modules. The data were retrieved from the 
Spitzer Science Center Data Archive using Leopard. The Spitzer 
IRS Custom Extractor (SPICE) was used to do the extraction of the 
spectra for each nod position from the 2D images. All spectra were 
cleaned for bad data points, spurious jumps and glitches, and then 
they were combined and merged to create a single 5-37yum spec- 
trum. We did not correct for flux offsets between different wave- 
lengths ranges but these are minor. Offsets can be caused by the 
fact that the different wavelength ranges use different slits, and our 
objects have sizes similar to or larger than the slit widths. 

Twenty-two of the objects were included in an HST SNAP- 
shot survey of 37 Galactic Bulge PNe (proposal 9356, PI Zijlstra). 
The WFPC2 images were taken in three different filters. Ho-, V, 
and [OIII], with typical exposure times of 2 x 100 sec, 60 sec and 
2 X 80 sec, respectively. The objects were placed on the CCD of 
the Planetary Camera (PC) with a pixel scale of 0'.'0455. Pipe-line 
reduced images were retrieved. Only the Hcf images are used here. 

The target selection for the original HST observation was 
based on a catalogued diameter less than 5". It included a few ob- 



jects without a known size. Other than this, the target selection was 
random among the known PNe in the direction of the Bulge. The 
fact that so many of the Spitzer sample were included in the HST 
observations indicate that the Spitzer target selection indirectly also 
made use of the diameter. The Spitzer samples were mostly selected 
based on 12yum flux. The overlap in selection reflects the brighter 
infrared emission from compact nebulae. 

VLT UVES echelle spectra were obtained for these 22 objects, 
using an exposure time of 600 sec. The spectra cover the wave- 
length range from 3300 to 6600A, with a resolution of 80,000. 
These spectra were used to investigate the nature of the central 
stars, based on the presence of stellar emission lines. 



3 RESULTS 

3.1 External causes: ISM interaction 

As argued above, PAHs could be present in an oxygen-rich en- 
vironment if they are acquired from the ISM. This can be done 
by gathering interstellar PAHs through the "wall" where the ex- 
panding nebula colli des with the ISM, creating a stationary shock 
jWareing et a"i]|2007h . in which the interaction region is dominated 
by ISM material. As the Bulge PNe have a high velocity dispersion 
with respect to the ISM, the "wall" is expected to be close to the 
nebula, and one-sided. 

To test this model, we checked the spatial location of the PAH 
emission bands, at 6.2 and 11.3yum, and compared these to that of 
the 12.8yum [Nell] line. The 12.8yum [Nell] emission line comes 
from the ionized region, so the position of the line traces the posi- 
tion of the PN and is not affected by extinction. On the other hand, 
if the PAHs have been swept up from the ISM, then their physical 
position would not be at the same position as the PN (or the [Ne II] 
for this case) so there should be an olfset in their positions. This 
was done for 40 objects in the sample. The peak of the [Ne II] line 
and the 11.3yum PAH feature were measured using the SH part of 
the spectrum. 

In a few cases we found evidence for an offset. In PN Hb 4, the 
faint 11.3yum feature is offset by about 6" from the [Ne II] 12.80yum 
emission. In M2-14, the faint 1 1.3yum emission is offset by 4" from 
the [Ne II] line. We also found an offset of 2" between the [Nell] 
line and the PAH bands at 11.3/zm in the case of Hl-16, Hl-50, 
H2-20, Ml-40 and Cnl-5. These smaller offsets may not be real, as 
the resolution at the position of the peak of the line and the feature 
is about 2". Hb4 provides the strongest case for a spatial offset. 
However, as shown below, it has at best weak evidence for mixed 
chemistry. 

Thus, there are only a few potential cases where interaction 
with the ISM, or other forms of ISM confusion, may contribute 
to the PAH features. In the large majority of cases, the dust and 
PAH features are co-located and centered on the nebula. In these 
cases the mixed chemistry is likely to have an internal origin, re- 
lated solely to the star and its ejecta. 

3.2 Internal causes: star and nebula 

3.2.1 Spitzer spectra 

Table[T]lists the 40 PNe we analyzed. We identified both PAH and 
silicate bands, as listed in the Table. PAHs bands were seen in all 
objects, except for 5 objects where only very faint 1 1.3/im/7.7yum 
bands and no other features were seen. Silicate emission was seen 
in 34 PNe. We found evidence for mixed chemistry (both PAH and 
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Table 1. Galactic Bulge PNe. Column 3 and 4 list the identified bands in the IRS spectra. Col. 5 lists the peak line flux of the 7.7/jm PAH band. Col. 6 indicate 
whether the star has emission lines [WC] or is of weak emission-line (wels) type. The last column indicates whether an HST image exists. A star in Col. 3 
indicates that the PAH bands are very faint and may be in some doubt. Mixed-chemistry objects are those unstarred objects which do not have "none" in col. 
4 (30 in total). 



PN G Name 


Name 


PAHsOjm) 


Silicates (fim) 


7.7/jm (mJy) 


Central star 


HST image 


OOS 3-01 1 


M 1-40 


6.2, 7.7, 11.3 


27.5, 33.8 


400 


wels 




002 2-09 4 


Cn 1-5 


6.2 7.7 11.3 


27.5 33.8 


280 






356 S+03 3 


Th 3-12 


6 2 7 7 8 6 113 


33.8 


180 






006.5-03.1 


H 1-61 


6.2, 7.7, 11.3 


27.5, 33.8 


150 


wels 


no^ 


006.4+02.0 


M 1-31 


6.2, 7.7, 11.3 


10, 27.5, 33.8 


140 


wels 




3^7 1-04 7 


H 1-43 


6 2 7 7 8 6 113 

U.^, 1,1^ O.U^ Xl.J 


10 27 5 33 8 


130 


[WCll] 


yes 


007.2+01.8 


Hb 6 


6.2, 7.7, 11.3 


27.5, 33.8 


120 




no 


358 7+05 2 


M 3-40 


6 2 7 7 8 6 11 3 


27.5 33.8 


111 






356 5-02 3 


M 1-27 


6.2, 7.7, 11.3 


27.5 


100 


[WCll] 


no^ 


004 Q+04 9 


M 1-25 


6.2 7.7 11.3 


27.5 33.8 


90 


[WC6] 




359 3-01 8 


M 3-44 


6.2, 7.7, 11.3 


27 5 33 8 


80 


[WCll] 




■XS'i Q+03 4 


H 1-19 


6 2 7 7 8 6 113 


27 5 33 8 


80 






000.1+04.3 


H 1-16 


6.2, 7.7, 11.3 


27.5, 33.8 


80 




no 




H 1-8 


6 2 7 7 8 6 113 


27 5 33 8 


70 






003.1+03.4 


Hen 2-263 


6.2, 7.7, 8.6, 11.3 


23 5 27 5 33 8 


60 




yes 


354 5+03 3 


Th 3-4 


6 2 7 7 8 6 113 

\J.^, I.I, O.U, 11. J 


10 27 5 33 8 


60 






359 9-04 5 


M 2-27 


6.2, 7.7, 11.3 


27 5 33 8 


60 


wels 


n^ 


000 0-06 8 


Hen 2-367 


6.2 7.7 11.3 


27.5 33.8 


60 






351.9-01.9 


K 5-4 


6.2, 7.7, 8.6, 11.3 


none 


50 




yes 


^5fS 9+04 4 


M 3-38 


6.2, 7.7, 11.3 


27.5, 33.8 


50 




yes 


002.8-01.7 


H 2-20 


6.2, 7.7, 11.3 


27.5, 33.8 


40 




yes 


006.8+04.1 


M 3-15 


6.2, 7.7, 11.3 


27.5, 33.8 


20 


[WC5] 




007.5+04.3 


Th 4-1 


7.7, faint 11.3 


10, 27.5, 33.8 


20 




yes 


^58 5-04 2 


H 1-46 


7.7, 11.3 


10, 27.5, 33.8 


18 




yes 


004.0-03.0 


M2-29 


7.7, 8.6, 11.3 


none 


10 




no 


355.9+03.6 


H 1-9 


7.7, 8.6 


none 


10 




no 


357.2+02.0 


H2-13 


faint 7.7, 11.3 


33.8 


8 




yes 


001.2+02.1 


Hen 2-262 


faint 6.2, 7.7, 8.6, 11.3 


33.8 


8 




no 


004.1-03.8 


KFL 11 


faint 7.7, 8.6, 11.3 


33.8 


5 




yes 


003.1+02.9 


Hb4 


faint 7.7 and 11.3 (*) 


27.5, 33.8 


5 


[WC4] 


no 


002.9-03.9 


H2-39 


faint 7.7 (*) 


23.5, 27.5, 33.8 


1 




yes 


006.3+04.4 


H2-18 


11.3 


33.8 







yes 


003.6+03.1 


M2-14 


faint 11.3 (*) 


27.5, 33.8 





wels 


yes 


356.5-03.6 


H2-27 


faint 11.3 (*) 


none 





[WCll] 


yes 


358.5+02.9 


Al 2-F 


faint 11.3 (*) 


33.8 







yes 


008.6-02.6 


MaC 1-11 


faint 8.6, 11.3 


33.8 







yes 


006.1+08.3 


M 1-20 


6.2, 8.6, 11.3 


none 





[WC8] 


no 


005.9-02.6 


MaC 1-10 


11.3 


33.8 







no 


008.2+06.8 


Hen 2-260 


8.6 


none 







no 


359.7-02.6 


H 1-40 


6.2, 11.3 


10, 27.5, 33.8 







no 



silicate emission) in 30 PNe, where we exclude those with dubious 
PAH bands. 

Table [T] is in order of PAH strength, where the PAH strength 
is defined as the peak line flux of the 7.7yum PAH band. From H2- 
18 onward, no 7.7^m emission feature is detected. The spectra are 
shown in the left-most panels of Figs.[T]|2l[3]and|4] We only show 
the short wavelength region (SL and HL) to cover the lOyum silicate 
feature and the PAH bands. The vertical red lines indicate the wave- 
length of PAH bands: 6.2, 7.7, 8.6 and 11.3yum . Offsets between the 
different spectral segments are visible in three cases. 

Th 3-4 (PN G354.5+03.3), H 1-46 (PN G358.5-04.2), and Th 
4-1 (PN G 007.5+04.3) show a strong broad band at lOyum, inter- 
preted as silicate emission. In all three cases, the peak is shifted 
to longer wavelengths, beyond lOyum. In the classification scheme 
of Sloan & Price (1998), they are of type SE4/5. The 13yum feature 
typically associated with the re ddest silicates (SEl-4) i n AGB stars, 
is not present in our objects. IPerea-Calderon et"^ show 



three more objects with silicate emission without PAH features - 
these also have reddish silicates, and are of type SE5-6, without a 
13yum feature. The reddening of the silicate emission may be at- 
tributed to large dust grains. 



3.2.2 Emission-line stars 

The right-most panels of these figures show the UVES spectra lo- 
cated around the strong stellar emission C IV-5800/581 lA located. 
This line is an indicator for [WC]-type stars. The y-axis scale is 
the same for all panels, except for two objects with strong emission 
lines. In the case of Hl-43 (PN G357. 1-04.7) we show the stellar 
broad line C III at 5696A together with the nebular narrow line 
of [Nil] at 5755A. T he UVES spectra confirm the classification 
of iGomv et id] ( |2004 but do not reveal any further emission-line 
stars. 
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Table [T] shows that the mixed chemistry phenom enon is not 
hmited to emission-line objects (cf. iGomv et alj2oToh . In our sam- 
ple of 30 PNe, 10 have emission-line stars and 20 do not. In the 
rejected sample of 10 PNe, 4 have emission-line stars and 6 do 
not. This confirms that in the Bulge, there is no clear relation be- 
tween mixed chemistry and class of central star. However, among 
the mixed chemistry objects, there is a strong tendency for the emis- 
sion line stars to have the strongest 7.7/jm band. Thus, although the 
mixed chemistry is itself independent of the stellar class, it is am- 
plified by emission-line stars. 

3.2.3 Morphology 

We use HST images of 22 of our objects to study their morphol- 
ogy. The Ha HST images are shown in the middle panels of Figs. 
[T] |2] [3] and |4l again in the order listed above (of decreasing PAH 
strength). The arrows indicate the orientation of the image and the 
bar indicates the spatial scale. The images are displayed using an 
"asinh" scale to enhance the contrast and to better see the faintest 
structures. (Absolute intensities are not used in this paper.) Note 
that the images of different objects are not to the same spatial or 
intensity scale. 

A few further i mages are available i n the literature. Ml -40 has 
a diameter of 20" jSchwarz et al.llT992h and is ther efore likely a 
foreg round nebula. Cnl-5 has a similar-sized halo dCorradi et al.l 
2003h and also is probably not in the Bulge. Ml -25 was imaged by 
Aaauist & Kwo3 ^19901) who find a torus with diameter 1.7". Hb 
6 and Hb 4 are als o not expected to be in the Bulge as their r adio 
fluxes are too high feiilstraet allll989l;lAaauist & Kwoldll990l) . 

The most important aspect of Figs.[Tl|2][3]and|4]is the pre- 
dominance of a very specific morphology: a dense torus with bipo- 
lar or multipolar flows. This appears to correlate well with the PAH 
band strength. In a few cases the torus appears to be seen pole- 
on (PN G358.7+05.2 and PN G002.8-F01.7), but in other cases the 
structure is easily visible. The confining tori have a range of sizes 
but are always well defined. The objects without 1 .Ijim feature, 
towards the end of the sequence, tend to show less structured mor- 
phologies. The clearest exceptions to the rule are PN G003. 6-1-03.1 
which has no 1.1 detection (and has a relatively large torus) and 
PN G003. 1+03.4 which does not show elongated lobes, possibly 
because of viewing angle. We note that one of the three silicate 
objects (PN G007.5+04.3) is unresolved by HST and is extremely 
compact. 

This type of morphology is not the most common in the 
Galaxy - most PNe are round or elliptical. In the lAC morphologi- 
cal catalogue (Manchado et al. 1996), only 10 bipolar/quadrupolar 
nebulae are present among 240 objects. Corradi & Schw arz (1995) 
find that 15% of Galactic PNe are bipolar. IZiilstral(l2007b finds 40% 
of compact Bulge PNe are bipolar or related. This is in contrast to 
our 70%. 

We conclude that the mixed chemistry phenomenon is 
strongly related to morphology, or the presence of a torus. There 
is no clear relation to the central star, apart from a possible ampli- 
fication of the PAH bands by emission-line stars. 

3.3 Modelling the Torus 

To test if the dust spectra can be related to the presence of a torus, 
we model the Spitzer spectrum for the representative PN Ml-31 
(FN G006.4+02.0 ), using the MC3D (MonteCarlo 3-D) program 
i Wolfetal.|[T999h . We did not incorporate PAHs in this model but 



Table 2. Model parameters for the PN Ml-31 



i?i„[AU] Rout[AU] Md„„(Mo) a /} ho[A\]] 
5000 50000 5.5x10-3 2.5 2 500 



only used astronomical silicate as the dust component as our aim 
is to reproduce the shape of the continuum due to emission from 
these dust grains. 

The model has 6 main parameters that can be adjusted to arrive 
at the best fit. Other parameters were kept fixed: we used an effec- 
tive temperature of the s tar of 57,000 K, a luminosity of 6,500Lo 
jGesicki & Ziilstrall2007l) . an inclination of the disk of 90 degrees 
(edge-on), with the distance to the Galactic Bulge taken as 8.5 kpc 
JSchneider & BucklevllT99^ . We did test calculations for a range 
of inclinations, but found that this parameter has a limited effect on 
the result, because of the limited opacity at lOyum of the torus. 

The Monte Carlo simulation calculates the radiative transfer 
for 100 wavelength points between 10"' and lO'^m, to calculate 
the dust temperature. It then uses 29 wavelengths between 1 and 
53/jm to calculate the SED. 

We obtained a best fit for the parameters listed in Table |2] The 
main parameters are: the inner radius of the torus (i?,,,), the outer 
radius (Rom), the dust mass (Mrf,„,), the exponent of the power law 
describing the density in the midplane (r ~°), flaring ( r^) and the 
scale height at 100 AU from /?„„ (ho) (Wood et aljl998l) . The fitted 
inner radius of the torus, of 5000 AU, is in good agreement with 
that measured from the HST image (0.55" which at 8.5 kpc is 4600 
AU). In Fig[5] we plot the Spitzer spectrum and the MC3D model. 
A good fit is obtained, with the model reproducing the steep onset 
of the continuum at 14jim. The model is a little below the measured 
flux at shorter wavelengths, which may be due to a contribution 
from PAH or very small grain continuum emission. We also calcu- 
lated a model with a mixture of silicate and carbon dust, but it did 
not reproduce the longer wavelength continuum as well, predicting 
too much flux longward of 24yum. 

The inner radius fits the HST image well but is relatively far 
from the star. This indicates that the fitted component is the ex- 
panding torus, and not a Keplerian disc. The model scale height is 
derived very close to the inner edge of the torus (100 AU outside 
of the inner radius of 5000 AU) where it is 10% of the radius, i.e. 
a thin torus. The flaring increases as r- and at the outer edge of the 
model, the height is equal to the radius. The HST image suggests 
a thickness of the torus which is intermediate between these val- 
ues. The dust mass indicates a gas mass of around 0. 1 M© (for an 
assumed gas-to-du st ratio of 200 ) The ionized mass has been esti- 
mated as 0.29 Mq jGesicki & Zii lstra 2007). This suggests that the 
modeled torus accounts for a large fraction of the nebular mass. 

Our calculation shows that the dust continuum emission can 
be explained as arising from the dense torus using oxygen-rich dust 
(silicate) to explain the SED: there is no significant carbonaceous 
component. The fitted torus is large and massive. This is a signifi- 
cant finding, as the traditional model for mixed chemistry poses a 
long-liv ed, compact oxygen -rich disc. Such a disc is in Keplerian 
rotation. IWatersetal.l(ll998h obse rved in the Red Rectan gle a low- 
mass disc of 500-2000 AU in size. lChesneau et d H2007l) observed 
a disc size of 9-500 AU in Mz-3 with a dust mass of IxIO^^Mq. 
We do not find evidence for such a component. 

In the traditional stellar evolution model, the silicates were 
ejected > lO'yr ago, with the star subseque ntly becoming carbon- 
rich prior to the ejection of the main nebula (lHabindll996[) . In con- 
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PN G 006.4+02.0 




PN G 357.1-04.7 



%0 
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57B0 SaOO 5B20 5340 
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uave length [A] 



5700 5720 5740 5760 
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G7B0 SBOO 5820 SS40 



Figure 1. CoiTelation image showing in the first column the IR Spitzer spectrum, showing the short wavelength region, the vertical red lines show the PAH 
bands. In the second column is the con'esponding HST image and in the third column we show a part of the UVES spectrum. 
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Figure 5. The Spitzer spectrum of the PN Ml-31. The dashed black Hne is 
the model calculated using the MC3D code. 



trast, the torus of Ml-31 is part of the main nebula and ejected 
< lO'' yr ago: at this time the ejecta were oxygen-rich. The mixed 
chemistry phenomenon is therefore unlikely to be related to the star 
becoming carbon-rich, as the chances of this occuiTing over such a 
short time scale are low, whereas the occurrence of mixed chem- 
istry appears to be relatively common. 

Interestingly, the two brightest PAH emitters in the mixed- 
chemistry sample, Ml-40 and Cnl-5, are foreground nebulae. 
These are emission-line stars for which the traditional model of 
a long-lived disc may hold, similar to objects such as IRAS 07027- 
7934 and the Red Rectangle. The mixed chemistry nebulae in the 
Bulge are a population which is distinct from these disc objects, 
with a different origin of the mixed chemistry. 



4 ORIGIN OF THE PAH IN BULGE PNE 

The previous analysis shows that the mixed-chemistry Bulge neb- 
ulae are likely oxygen-rich. This is consistent with expectations 
based on the stellar population in the Bulge. The Bulge is be- 
lieved to cont ain a lOGyr old pop ulation, with little or no trace of 
younger stars jZoccali et al.l2003l) . although recently, some indica- 
tions have been fou nd for a younger, highly metal-rich population 
teensbv et allboiol) . 

Carbon stars originate from third dredge -up, which requires 
a mi nimum stellar mass of M, >~ 1.5 Mg jVassiliadis & Woo3 
1 1993 ). It is not expected to occur in lOGyr old populations. Fur- 
thermore, at high metallicity the formation of carbon stars is in- 
hibited by the large oxygen abundance which needs to be over- 
come b y the primary carbon. Carbon stars are rare in the inner 
Galaxy i Le Bertre et al.l I2OO3I) and absent in the Galactic Bulge 
( lFeastll2007hT consistent with this expectation. 

This strengthens the case that the mixed chemistry is occurring 
in oxygen-rich gas, expelled by oxygen-rich stars. 

We find a strong correlation between mixed chemistry and 
morphology, specifically, the presence of a massive torus. This in- 
dicates that the chemical pathway to the PAH molecules occurs 
within their high density regions. These tori remain molecular for 
some time after the ionization of the nebula has started, due to 
the trapping of the ionization fr ont and due to the shielding ef- 
fects of the dust column density jWoods et al]|2005l) . A mini PDR 



(photon-dominated region) is expected in these molecular regions 
jPhillips et ai]|20IOl) . slowly being overrun by the advancing ion- 
ization front. 

A mini-PDR within a dense torus is a plausible location for the 
PAH formation. In interste llar clouds, PAH emission is often as so- 
ciated with PDR regions ( Joblin et alJbOldlKassis et alj|2006l) . It 
is unclear whether PAHs form in these regions or whether they be- 
come visible due to UV excitation. In our objects, because of the 
close relationship between morphology and PAHs, it seems likely 
they form in the regions where they are observed. 

In AGB stars, the dichotomy between O-rich and C-rich chem- 
istry can be broken, as is shown by the observation of hot water 
in the carbon stai-s IRC -I-I02I6 jPecin et alj|2010l : iNeufeld et all 
l20IlHCherchnefj|20Ilh . This may be due to collisional destruc- 
tion of CO in the shocks 1 Duari et al .I1999I) or interstellar UV pho- 
tons d Agiindez etai]|201oir The latter authors find that for O-rich 
circumstellar envelopes, the penetration of interstellar UV photons 
into the inner layers leads to the formation of molecules such as 
CH4 and HCN with abundances relative to H2 in the range of 10"* 
to 10"^. This chemistry is initiated by CO dissociation; the dissoci- 
ation energy of CO is 1 1 . 1 eV. 



4.1 Chemical reactions 

Here we investigate whether a C-rich chemistry can be similarly 
driven by UV photons penetrating the dense, oxygen-rich torus. 
The penetration is limited by the dust extinction: in our HST image, 
internal extinction is clearly visible in M3-38 (PN G356. 9+04.4) 
which has the most compact torus. 

Our model, which contains gas-phase chemistry only, is based 
on that of Ni Chuimin (2009) (see also Millar et al. 2007), devel- 
oped for conditions in PDR regions. I t uses an extrapol ation of the 
Meudon 2006 PDR chemistry code I Le Petit et alj[2006) , and as- 
sumes a semi-infinite and constant-density slab with a density of 
2 xlO"* cm"', UV radiation enhanced by a factor of 60 over the 
average interstellar field, treats thermal balance self-consistently, 
and contains a mixture of reactions ext rapolated from the M eudon 
che rnical files and also th e UDFA 2006 jWoodall et alj2007lfl Rate 
99 ( lLeTeufreta"i]|2000h and Ohio State OSU-O^EI databases. In 
Fig. |6]we show the behaviour of the temperature against A, . As 
expected, the temperature does not vary much between the three 
models with each displaying an increase in cooling by atomic O 
as its abundance increases. Formation and destruction of large hy- 
drocarbon chains, that is with more than 4 carbon atoms, of the 
form C„H,„, (n < 23, m < 2) and their associated ions were added. 
This resulted in a chemical network of 269 species connected by 
about 3500 reactions involving ion-neutral and neutral-neutral col- 
lisions, photodissociation and photoionization, dissociative and ra- 
diative recombination with electrons and radiative association reac- 
tions. The goal of the models is to study the abundance s of large hy- 
drocarbons. The specific formation of benzene rings (IWoods et al.l 
I2002ll2003h and PAHs is not included, but cyclisation is assumed to 
be possible for chains of 23 carbon atoms. Laboratory spectroscopy 
has shown that this does not yet happen spontaneously for chains 
of 13 carbon atoms ( lGiesenlll994l) . 

The elemental carbon abundance is kept constant at the stan- 
dard value of 1.32x10"* relative to the total hydrogen abundance 
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Figure 6. Temperature as a function of visual extinction in our three models. 
The temperature is calculated self-consistently in the constant-density slab. 



Table 3. Example of model reactions of different types. 
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Figure 7. Calculated abundances for O, C, C-f and CO for the three models 
with varying C/O ratios. Only the O atom abundance varies, other abun- 
dances are almost identical in the models and are shown here for the case 
of C/O = 1/2. 
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dSavage & SembacHI 19961 ) while the elemental oxygen abundance 
is varied to give a range of C/O values of 2/3, 1/2 and 1/3. 

The formation pathways for the carbon chain species in the 
chemical model include neutral-neutral C atom reactions, radia- 
tive association reactions with carbon atoms, and ion-molecule re- 
actions with C^. Destruction of these chains occur in reactions 
with oxygen atoms (O + C„), photodissociation and photoioniza- 
tion, dissociative recombination and electronic recombination reac- 
tions (Table [3](. The oxygen-destruction reaction rat es are adopted 
and e xtrapolated from the Meudon 2006 chemistry dLe Petit et al.l 
l2006h . The oxygen destruction rates for Ct and C4 are assumed 
to be reasonably fast (1.5 x 10"" and lO^'^cm'.?"', respectively), 
while the reaction wit h C3 is assumed to have an energy barrier 
( IWoon&Herbstlll996l) . 

Example reactions of various types included in the model are 
listed in Table[3] The last column shows the source of the adopted 
rate coefficients, either from the Meudon model, the UDFA, or 
OSU-03 databases. 



4.2 Results 

If PAH emission does originate from the torus and not swept-up in- 
terstellar gas as is indicated by our observations, then it is valid 
to consider whether the PAH themselves can be formed in this 
oxygen-rich, UV-irradiated environment. Since there is no detailed 
chemical kinetic model for PAH formation, we may use our hydro- 
carbon chemistry as a proxy for PAH formation, noting that even 
the presence of abundant hydrocarbons in a PDR would not neces- 
sarily prove that PAHs are formed in the same environment. 

The PDR model calculations show that the chemistry first 
tends towards complexity around the C+/C/CO transition at Ay ~ 
1.4 mag. This transition is seen in Fig.|7]which plots the abundances 
of O, C, C-l- and CO as a function of A, . As expected, only O shows 



a difference between the three models, for the other species we sim- 
ply plot the results for the model with C/O = 1/2. 

Figs. [8] |9] and [TO] show the calculated fractional abundances 
(relative to H2) for the molecules C2H, CgH and C23H2, respec- 
tively, as function of A,,, which measures depth into the slab, and 
C/O ratio. 

Two hydrocarbon peaks are seen, one at Ay ~ 1.5, and one 
at Av ~ 4-6. The first peak occurs in the region in which C^ and 
H2 are abundant with the hydrocarbon chemistry initiated by the 
C^-H2 radiative association reaction. Subsequent reactions with C 
and C^ form CiH. C2H is the main seed molecule for the forma- 
tion of large r hydrocarbons, as it is in the carbon-rich envelope of 
IRC+ 10216 jMillar et al.l200(]h , and reactions continue to build up 
to CgH, although large(r) hydrocarbons are not highly abundant at 
low extinction. 

The second peak, at higher Ay ~ 4 - 6, occurs where atomic 
carbon becomes abundant and drives chemistry through proton 
transfer with Hj to form CH+ which subsequently reacts with H2 to 
form simple hydrides. Here the chemistry produces a surprisingly 
large abundance of the biggest hydrocarbons in the model. 

Oxygen atom destruction reactions are the dominant loss of 
larger hydrocarbons at high Ay where there is little C*. The re- 
action rates of C2H and C with O are rather uncertain: they dif- 
fer by an order of magnitude between the UDFA and the Meudon 
chemical models. However, because of the importance of O, the 
abundances depend on the C/O ratios: models with less O produce 
higher abundances of large hydrocarbons. The difference amplifies 
for the larger molecules: whereas for C2H and CgH the difference in 
abundance between C/0=l/3 and 2/3 is within an order of magni- 
tude, for C23H2 it is four orders of magnitude. Even a small increase 
in the C/O ratio makes a large difference here. 



4.3 Discussion 

Under the assumptions of the chemical model, it is clear that larger 
hydrocarbons can form within oxygen-rich environments. It is as- 
sumed here that a mechanism exists to convert these long linear 
hydrocarbon chains into PAH molecules. Such a discovery is sur- 
prising, but clues have been previously found in the inner regions 
of oxygen-rich stars, where small carbon-bearing molecules have 
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Figure 8. Calculated abundances for varying C/O ratios for C2H. 
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Figure 9. Calculated abundances for varying C/O ratios for CgH. 



been observed, and explained by means of shock chemistry (Cher- 
chneff 2006) or penetration of UV photons (Agundez et al. 2010). 
Furthermore, in the oxygen-rich environment of a protoplanetary 
disc. Woods & Willacy (2007) showed that benzene and large 
polyynes of up to C8H2 can form through ion-molecule chemistry 
in suitably dense conditions, where cosmic- and X-rays provide a 
source of ionisation. In this latter model, and in the ISM (McEwan 
et al. 1999), and in proto-PNe with carbon-rich tori (Woods et al. 




Figure 10. Calculated abundances for varying C/O ratios for C23H2. 



2002, 2003) closure of the aromatic ring occurs through the reac- 
tion of small (Cj 4113^5) hydrocarbons, which are not included in 
our model due to the simple nature of the Meudon chemistry (Table 
O. However, these small hydrocarbons are typically formed by the 
destruction of larger polyynes, and through reaction with atomic 
carbon, which are abundant in the chemical model. Once this first 
aromatic ring closure has occurred, subsequent aromatic growth is 
rapid in comparison (Cherchneff et al. 1992, Wang & Frenklach 
1997). 

If we assume that the underlying abundance distribution of 
hydrocarbon chains is related to that of PAHs, then we can make 
some inferences on the properties of PAH emission in the Galactic 
Bulge PNe. First, a reasonable extinction is needed. At low Ay (<2), 
where carbon is ionized, the abundances of the longest chains are 
strongly suppressed by photodissociation and photoionization. This 
still allows for the formation of shorter linear hydrocarbons, but 
PAH formation is not expected. Efficient PAH formation requires 
sufficient extinction to render carbon neutral and to minimize the 
effects of photodestruction. The optimum extinction in our model 
is around Ai/ ~ 4. At even higher extinction, carbon is locked up 
in CO and fewer of the longest hydrocarbon chains form, although 
intermediate chains (CgH) still form. 

In order to generate sufficient reactive C and UV photons 
are required to dissociate CO contained in the dense molecular 
torus. Thus, the central star needs to be hot enough to generate 
these photons. Our models use a radiation field of 60 times the in- 
testellar radiation field which corresponds to a star of spectral type 
B, or temperature around 15 OOOK. This suggests that the PAHs will 
not form until the star is hot enough to begin to ionize the nebula. 
By this time, the nebula will have expanded and the extinction has 
dropped considerably. 

The need for the torus may therefore be to provide an irradi- 
ated region dense enough to allow a photon-dominated chemistry, 
in a nebula old enough to allow for an ionizing star. Such a require- 
ment would explain the close relation between PAH emission in 
oxygen-rich PNe, and the presence of a dense torus. This same re- 
lation also holds in non-Bulge PNe, e.g. NGC 6302 and Roberts 22 
( IZiilstraetal]|200lh . 

The predicted abundances show a strong dependence on 
C /O ratio (Figs. [8]_|^ HOt . This may explain the finding of 
ICohen & Barlow! ( |2005|) that PAH emission at lower C/O ratio is 
faint. An increase in the C/O ratio, even if it remains below unity, 
would benefit PAH formation. But such an increase is not expected 
in Bulge PNe which do not experience carbon dredge-up. 

Finally, even though this model goes some way to explaining 
the observed mixed chemistry, it remains unexplained why PNe 
with [WC] central stars would have stronger PAH emission. Part of 
this correlation may be due to the fact that the brightest objects in 
Table[T]contain some foreground nebulae. However, even excluding 
these, it still appears that the emission-line stars dominate the top 
end of the Table. 

Two possible explanations can be envisaged. First, the stellar 
winds may inject some carbon into the surrounding nebulae. This 
would locally enhance the C/O ratio, and increase the PAH foma- 
tion rate. 

The second explanation could be r elated to the energy d eposi- 
tion by the stellar wind into the nebula. iGesicki et aO ( |2006|) show 
that nebulae around [WC] stars are turbulent indicating that some 
energy deposition into the nebulae occurs, for instance from the 
stellar wind. The turbulence may aid in dissociating CO, or more 
likely, causes compressed regions with a larger range of extinctions. 
The second explanation appears more plausible. However, the rela- 
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tion between stellar winds and PAH emission needs to be studied 
further. 

We finally note that the Ay in Figs. [H H] [TO] [7] is time- 
dependent. As the torus expands, the extinction drops. The local 
conditions thus evolve from right to left in the Figures. This pre- 
dicts two layers within the torus, one with the smaller hydrocar- 
bons at Ay ~ 1.5 and one with PAHs at Ay ~ 4, which are both 
moving inward into the torus as the nebula evolves. This makes the 
layer at lower extinction, with shorter hydrocarbons, a product of 
the destruction of the PAHs by the increasing UV radiation field. If 
the PAHs are more UV-stable than the largest chains in the model, 
then PAHs could be present at lower extinction than predicted here. 

5 CONCLUSIONS 

We summarize the main conclusions of this research in the follow- 
ing list: 

• We analyzed observations of 40 PNe towards the Galactic 
Bulge using the IRS instrument on board Spitzer. We found that 
30 of them present a mixed chemistry phenomenon. From these 40 
we used HST observations of 22 to analyze their morphology. 

• A strong correlation was found between strength of the PAH 
bands and morphology, in particular, the presence of a massive 
torus. Data from VLT/UVES was also used to study their central 
stars. We find that there is no correlation between the classifica- 
tion of the central star and the mixed chemistry phenomenon, al- 
though PNe with emission line stars (wels or [WC]) tend to have 
the stronger PAHs emission. This could be related to carbon injec- 
tion or energy deposition by the wind into the nebula, enhancing 
PAH formation. 

• We modelled the Spitzer spectrum for the representative PN 
Ml -31. The SED can be explained using oxygen-rich (silicate) 
dust. The fitted torus is large and massive, with a size of 5000- 
50000 AU and a mass of 5xlO"'Mo. Because of the age of the 
torus of Ml-31 (< 10''yr) the mixed chemistry phenomenon is un- 
likely to be related to a change of composition of the stellar ejecta 
over time. Such a change would be unexpected in the Bulge pop- 
ulation (which does not show third dredge-up), and would be very 
rare over such a short period in any case, whereas the occurrence 
of mixed chemistry appears to be relatively common. 

• For the chemical analysis of the PAH formation in an oxygen- 
rich environment, we used the gas-phase models from Ni Chimin 
(2009). These are an extrapolation of the Meudon 2006 PDR chem- 
istry code to include the chemistry of large hydrocarbons. They 
cover a range of A,, and C/O ratios of 1/3, 1/2 and 2/3. Although 
the model does not include PAH formation, it shows that large hy- 
drocarbons can form in oxygen-rich environments and that their 
abundances can be relatively large in moderately shielded regions. 
The model predicts two layers, one at Ay ~ 1.5 where small hy- 
drocarbons form from reactions with C^, and one at Ay ~ 4, where 
larger chains (and by implication, PAHs) form from reactions with 
neutral, atomic carbon. 

• The torus observed in all of the PNe that present the mixed 
chemistry phenomenon could provide the dense, irradiated envi- 
ronment needed to form the PAHs detected. 
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